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In the vicinity of the onset of an instability, we investigate the effect of colored multiplicative 
noise on the scaling of the moments of the unstable mode amplitude. We introduce a family of zero 
dimensional models for which we can calculate the exact value of the critical exponents /3 m for all the 
moments. The results are obtained through asymptotic expansions that use the distance to onset as 
a small parameter. The examined family displays a variety of behaviors of the critical exponents that 
includes anomalous exponents: exponents that differ from the deterministic (mean-field) prediction, 
and multiscaling: non- linear dependence of the exponents on the order of the moment. 



I. INTRODUCTION 

Critical exponents are usually introduced in the con- 
text of continuous phase transitions at equilibrium. The 
order parameter (for instance the magnetization for a sys- 
tem of spins or the density difference between the phases 
for the liquid-gas critical point) depends on the distance 
from the critical point as a power-law. The exponent of 
the power-law, traditionally named /?, is one of the criti- 
cal exponents of the system. Mean-field approximations 
simplify the analytical approach and allow to calculate 
(3. In this case simple rational values are obtained that 
depend on the nonlincarity (itself usually constrained by 
the symmetries) of the system. However, because of ther- 
mal fluctuations the mean-field results are not always 
correct for low spatial dimensions. In particular j3 can 
take non-mean-field values [l[ . Apart from the few cases 
in which exact analytical results exist, renormalization 
methods are used from which /3 is obtained as a scries in 
the critical dimension (the dimension above which mean- 
field results apply) minus the spatial dimension 0. The 
importance of spatial dimensionality appears also clearly 
in the fact that equilibrium phase transitions do not oc- 
cur when the spatial dimension is too small and short 
range interactions are considered. For continuous phase 
transitions at equilibrium, thermal fluctuations, nonlin- 
earities and spatial variations must be taken into account. 

In these systems at equilibrium, fluctuations are addi- 
tive terms in the equation for the order parameter. In 
out of equilibrium systems, fluctuations can be coupled 
differently to the order parameter and in several models, 
multiplicative noise is considered. Then, the analogous of 
phase transition can occur even when space is not taken 
into account and only the time evolution of a finite num- 
ber of modes is considered. Since spatial dimensions are 
not taken into account we refer to these models as zero 
dimensional. The simplicity of zero dimensional models 
allows a much more thorough analytical investigation and 
in a few cases the calculation of the full probability dis- 
tribution function (pdf) is possible. This is the direction 
that we pursue in this work. 

We introduce a family of zero dimensional models for 
which wc can calculate (in the small deviation from crit- 
icality limit) the stationary pdf of the system and thus 



we can obtain the exact value of the critical exponents 
Pm for all the moments. Despite the simplicity of the 
models the results show a rich behavior. The critical ex- 
ponents can differ from their deterministic values and in 
some cases vary continuously with the system parame- 
ters. In the later case a non-linear dependence of the 
exponents on the order of the moment is observed and 
thus the system displays multiscaling. 

In section [IIJ we present what is known in the deter- 
ministic limit and in the case of a white noise. The family 
of models under study is presented in section IIIII In sec- 
tionllVland sectionlVl we present two particular cases the 
second of which results in anomalous exponents. In sec- 
tion [VH our results arc interpreted based on a heuristic 
arguments from which the value of the critical exponents 
is understood. We conclude in the last section. In 0], 
we reported on the value of the exponent of the first mo- 
ment obtained for certain values of the parameters. The 
associated asymptotic expansion is presented in detail 
here together with several new expansions that are valid 
for other parameter values. Overall we are now able to 
calculate the whole set of exponents of all orders of the 
described model. 



II. ZERO DIMENSIONAL BIFURCATIONS 

We consider the evolution of an order parameter, x, 
which is a function only of time t. It satisfies the 
Langcvin equation 



X = /ax — \x\ n x + £x, 



(1) 



where fj, is the parameter that controls the instability, 
n > characterizes the nonlinearity (for instance n = 2 
for cubic nonlinear terms, n = 4 for quintic ones) and £ 
represents random fluctuations with zero mean (noise). 
From now on, in the case that £ is white noise, we use 
the Stratanovich interpretation [J] . We note that the so- 
lution x(t) conserves its sign and we thus restrict only to 
positive values for x. We are going to characterize the be- 
havior of x using its moments evaluated in the long time 
limit that we write here as (x m ), where (•) stands for av- 
erage over the realizations of the noise. The stability of 
the x = solution is determined by calculating the value 
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of the growth rate 7 = (x/x) for the linear system Q. 
Here 7 = /z thus the onset of the instability is given by 
fi = 0. For /i < the only attracting solution of the sys- 
tem is x = and thus all moments are zero. For positive 
values of /x, (x m ) takes non-zero values whose amplitude 
has a power-law dependence on /1: (x m ) oc yfi m . The ex- 
ponents of these power-laws f3 m are of primary interest 
in this work. Explicitly we define 



A, 



lmlog(<0)/log(/i). 



(2) 



In the deterministic limit (£ 



0) and for positive \i 
the long time solution satisfies limt_j. 00 x m = /x m / n . The 
critical exponents are thus (3 m = ^ . We will refer to this 
scaling as mean-field or deterministic scaling. 

A second well-studied limit is obtained when £ is a 
Gaussian delta-correlated noise, i.e. (£(t)£(t')) — 25(t — 
t'). Then the stationary probability distribution function 
(pdf) of x, P(x), satisfies the one dimensional Fokker- 
Planck Equation 



d x (fj,x 



P + d x xd x xP = 0. 



Its solution is given by 



P(x) 



n l-fi/n 

T (M/n) 



(3) 



(4) 



where the normalization condition J P(x) dx = 1 has 
been used. The moments can then be calculated as 



= / Px m dx 



(5) 



that, due to the singularity of P at x — 0, result to 
(x m ) ex /j, in the small jU-Hmit. We thus have (3 m = 1 for 
all moments m. This property is an effect of the noise on 
the dynamics of x close to the onset. Indeed, the time 
series alternates between phases where the value of x is 
cither large and nonlincarities are importantfon-phascs) 
or it fluctuates close to zero (off- phases) [a-Q. This be- 
havior is called on-off intermittency and an example of 
time series is displayed in fig. [TJ The mean duration of 
the off-phases, say T OFF can be estimated by considering 
the evolution of z = log(:c) for which eq. ([1} is written 
as 



(6) 



Thus z displays Brownian motion with a small drift (p) 
when x <C 1 while it is repelled towards the origin x = 
by the nonlinearity when x is order 1. As a result the 
duration of the off-phases diverges as T OFF ~ fiT 1 while 
the duration of the on-phases remains finite. During the 
on-phases, x achieves finite values, say x NL that do not 
depend on /1 (in the small /i limit). An estimate of the 
moments is given by (x m ) ~ T ON x™ L /(T ON +T OFF ) from 
which the linear dependence of the moments on fi is re- 
covered @. 



III. BIFURCATIONS IN THE PRESENCE OF 
COLORED NOISE 

For colored noise, it has been shown that the regime 
of on-off intermittency is controlled by the value of the 
noise spectrum at zero frequency, D = J Q (£(t)£(0))<fi 
0. As long as D is non zero, the behavior for very 
small /i is on-off intermittency [9(- In what follows, we 
examine the properties of the bifurcation when the noise 
has vanishing spectrum at zero frequency. 

We consider the family of models 



V 



fix - x n+1 
lt-Fy(y)} 



■X[Z-F V (V)] 



(7) 

where £ is a Gaussian white noise, (£(£)£(£')) = 2<5(i — 
if). We have introduced the potential F(y) which is a 
function of y only, and F y its first derivative. 

The amplitude of the order parameter, x, undergoes a 
bifurcation at fi = and is subject to a multiplicative 
noise y = £ — F y {y). Provided the stationary distribu- 
tion of y has a finite second moment, the spectrum of y 
vanishes at zero frequency. Indeed, for initial conditions 
2/(0) = 0, we can write 



dt 



2yy = 2y(t) / y(t')dt' . 
Jo 



(8) 



By averaging over the realizations and taking the long 
time limit, the last expression is the integral of the au- 
tocorrelation of y. This is also its spectrum at zero 
frequency using the Wiener-Kintchin theorem. At long 
time, if the second moment of y tends to a constant, then 
y has vanishing spectrum at zero frequency. 

The stationary pdf for y, Tl(y), satisfies the equation 
CqII = 0, where the linear operator Cq is defined as 



c Q u^d y [(d y F(y))n] + d^u, 



(9) 



and the normalization condition J Tidy = 1 is assumed. 
Equation CqII — has the solution 



U(y) = Exp[-F(y)]/N . 



(10) 



where N is a normalization constant. 

To isolate the noise term we make the following trans- 
formation w = log(x) — y. The Langevin equation be- 
comes 



w = /i 



n(w+y) 



(11) 

(12) 



The Fokker-Planck equation for the stationary joint 
pdf P(x, y) in these coordinates then reads 

a* Ox - e n{ - w+ ^)P = d y [F y P] + d 2 y P = CoP (13) 

Solving the partial differential equation (TT3|) for all values 
of the parameters is out of reach. Since we are interested 
in the critical behavior, we successively introduce various 
asymptotic approaches in order to determine the critical 
exponents. 
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(a) (b) (c) 

FIG. 1. Time series of the solution of eq. (0 for /i = 0.01. Top: linear scale and bottom in log scale, a) F — corresponding 
to a white noise; b) F = "fY 2 /2, Ornstein-Uhlenbeck noise with 7 = 1.5; c) F = v\Y\, with v = 0.75. Note the differences in 
the y-coordinate values. 



IV. STEEP POTENTIAL - RECOVERY OF THE 
DETERMINISTIC SCALING 

A. The general case 

The top panel of figure [5] shows the location in phase 
space of 10 3 trajectories for different values of [i (differ- 
ent colors) obtained by the numerical integration of the 
Langevin equations [7] for a steep potential that is de- 
fined in section IV. B. It can be seen that as \i becomes 
smaller the distribution is concentrated around a value 
of w = wo that depends on \i. Inspired from the nu- 
merical results and without any assumptions yet on the 
functional form of F, we make the following change of 
variables w — wq + ln(/x)/n + Uy/JL where wq = ln(Xo) 
is a constant that will be set by the expansion. The 
Fokker-Planck equation then reads 



C P = VMcyi - X%e<^ u+ ^)P. 



(14) 



We expand the p.d.f. as P = P + y^/ZPi + \iP-i + . . . . 

To lowest order we obtain the equation for the station- 
ary distribution in y 



CvPo = . 
The solution of which can be written as 
P = A(u)Tl (y) = A(w)e" J 



(15) 



(16) 



where the amplitude A(u) is left undetermined. To next 
order we have 



C0P1 = d u (l - X£e ny )P . 



(17) 



Integrating this equation over y provides us with a 
solvability condition. More precisely, since CqP\ is a to- 
tal gradient, integration over y makes the left hand side 
equal to zero and we are thus left with 







(1 - X^e ny )e~ F dy, 



This condition determines the value of Xq to be 

l/n 



X 



J^e^dy, 



(18) 



(19) 



We can then solve for Pi and we obtain 



Pi = A u e 
= A u U 1 (y), 



V" / {1-X^e nv ')e- F 'dy'dy" . 

Joe 



(20) 



where F' and F" denotes that the function F depends 
on the variables y' and y" respectively. To second order 
we have 

C P 2 = d u [(1 - X o n e"*0Pi - uXge^Po] . (21) 

Using again the solvability condition, we integrate over y 
and obtain 

/+00 r+ao 
(1 - Xge^ILidy - AuX™ / e^~ F dy = , 
-CO J —00 

(22) 

that leads to A — e _<5 " 2 with 



-X" 



+OO 



a ny-F(y) 



+00 



dyj I (X^e n y-l)IL 1 (y)dy. 

J — OO 

(23) 

Denoting Q{y) = j V _ 00 {X%e ny ' - l)e~ F ' dy' , we can show 
by integration by parts that 



+00 



{x^y - ljiir 0/) 



+00 



->2 F 



e 1 'dy>0. 



thus S > 0. The zeroth order solution then becomes 
1 



N 



(24) 



where N = ^pj e F dy. The moments can be cal- 
culated as 



N 



rrt/ft 

/ e ray-\-7n^/Jlu-8u 2 -F{y) fiyfa 

J —00 



.(25) 

The integral and N are independent to first order in fi. 
Provided the integral and X are finite, the moments 
follow the scaling (x m ) oc \i m l n . Therefore, we recover 
the deterministic exponents j3 m = m/n. 
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B. An example: y is the Ornstein-Uhlenbeck 
process 

A simple potential for which the former expansion is 
valid consists in F(y) = \^y 2 ■ This case corresponds 
to the Ornstein-Uhlenbeck process for the variable y [3]. 
Time series of a; are presented in fig. [TJ As mentioned, we 
show in fig. [2] the position of several trajectories in phase- 
space from which the concentration of the p.d.f around 
the value wq appears clearly. In this case the p.d.f. at 
lowest order becomes 



Pn = 



1 



: exp 



& i 1 i r \\2 1 2 

— {w-wo m(/i)) --71/ 

/i n Z 



, (26) 



with wq = logAo = —71/(27) an d 5 is given by the inte- 
gral (|23|) . Note that Xq takes very small values when 
7 is small which requires high accuracy when solving 
numerically the Langevin equation ([7]>. In fig. [5J wc 
present the marginal probability J P{uj 1 y)dy as a func- 
tion of (w — wq)/^ 1 / 2 for different values of /z and 7=1. 
The numerically computed pdf agree well with the theo- 
retical expression [2U 



C. Crossover: Two non-commutative limits 

Setting 7 — in cq. ([7]), we recover eq. <jTJ) with a white 
noise and this is the case in which on-off intermittency 
takes place. In this case, the exponents are /3 m = 1. This 
is at odd with the deterministic scaling (3 m = m/n that 
we have predicted for non-zero 7 (but possibly arbitrarily 
small) . 

Some insight on this problem of exchange of limits can 
be obtained by investigating the small 7 limit (or small 
amplitude of F in general). In this case, we write P — 
Po+7 1 ^ 2 Pi +7P2-K • ■ and use the fast z = w + y = ln(x) 
and slow Y = 7 1 / 2 ?/ variables. The Fokker-Planck in 
these coordinates reads 

c^(/j-e n2 -7 1/2 y)P = d Y iYP+d 2 P+jd Y P+2 1 1/2 d z d Y P ■ 
At lowest order in 7 we get 

9 z (/i- e " z )Po-^P = 0, (27) 

that leads to 

P = A(Y) cxp[/jz - e nz /n] = A(Y)U{z) . 

Note that II is the p.d.f. obtained when the noise is 
white, i.e. associated to on-off behavior. At next order 
we have 



3 2 (/z 



)P 1 -d 2 P 1 =d z (YP Q + 2d Y P 



(28) 



The solvability condition (here integration over z) does 
not set the amplitude A(Y) at this order, however eq. 
can be easily solved to obtain Pi . 
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FIG. 2. Top panel: Position in phase-space of 10 3 realiza- 
tions of the numerical solution of the Langevin Equation 
with F — ^"fy 2 , and 7 = 2. Different shades/colors (on- 
line) correspond to different values of as indicated. The 
dashed line corresponds to w + y — where the nonlinear 
term is of order one. Lower panel: The marginal probability 
n( w ) = / P{y> w )dy from the numerical investigation (dots) 
and the analytical prediction (solid line). Shades and colors 
are the same as in the top panel. 



At second order we have 

d z ( f i-e nz )P 2 -d 2 P 2 = d z (YP 1 +2d Y P 1 )+d Y YP +dyP . 

(29) 

Integrating over z, we obtain 

= d Y (AY + dyA), 

that leads to 

P =cxp[ i uz-e"7n- 72/ 2 /2] 



(30) 



(31) 



Using this result to calculate the moments for small /z, 
we recover the on-off scaling j3 m = 1 for all moments m. 
The validity of this expansion however holds only when 
7 tends to zero with fixed fi and thus it does not provide 
the actual critical exponent. The two limits of small (i 
and small 7 cannot be exchanged: 

hm/3 m ? lim hmlog(( a; m ))/log(/ / ). 
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FIG. 3. First moment as a function of fj, obtained from the 
numerically computed solutions of eq. [7] for F = 7j/ 2 /2 and 
five different values of 7. The first moment has been rescaled 
by the deterministic scaling fj, 1 ^ 2 . 



This explains the crossovers that are observed if we 
calculate the moments numerically for small values of 7. 
In fig. [3j we display the first moment as a function of fj, 
for n = 2. For large 7, the mean-field exponent (3 = 1/2 
is found for \i up to unity. In contrast for small 7, this 
exponent is recovered only for very small fx. For larger 
fi, an apparent exponent /3 smaller than unity can be 
estimated which traces back to the on-off behavior and 
is only a cross-over. We note that in an experiment or 
in a numerical simulations, such cross-overs can easily 
be interpreted as anomalous exponents. Indeed they can 
be observed on a large range of \x and the deterministic 
result is only recovered for very small \x. 



6 
A 

2 



-2 

-4 
-6 

-12-10 -8 -6 -4 
w 




-2 



1 .00 



0.10 



/x=5.60 10" 
© ^=3.16 10 
□ ,14=1.78 10 




0.0' 

-6 -4-2 2 

w— In(/i.) / v 

FIG. 4. Top panel: Position in phase-space of 10 3 realiza- 
tions of the numerical solution of the Langevin Equation [7J 
with F = v\y\, and v = 1. Different shades/colors (on- 
line) correspond to different values of /1 as indicated. The 
dashed line corresponds to w + y = where the nonlinear 
term is of order one. Lower panel: The marginal probability 
Ii(w) — J P(y,w)dy from the numerical investigation (dots) 
and the analytical prediction (solid line) 



V. ANOMALOUS EXPONENTS FOR A LESS 
STEEP POTENTIAL 

The expansion presented in section IV. A fails when the 
denominator in eq. (|19[) diverges. A simple potential F 
for which this expansion can break down is F = v\y\. 
In this case y follows a Brownian motion with solid fric- 
tion (l~5j . Time series of x, displayed in fig. [IJ show an 
intermediate behavior between on-off intcrmittency and 
the one described in section III.B. Phase-spaces for the 
numerically computed solutions of eq. ([7]) arc displayed 
in fig. gj 

For this choice of F, we need to consider two cases 
separately. For v < n, the expansion of section III breaks 
down at lowest order and a different asymptotic must 
be performed. For v > n, the expansion remains valid 
for the first moments but need to be modified for larger 
moments. 



A. The case v > n 

In this case the previous expansion works only for the 
calculation of the moments of order smaller than v. To fix 
this we need to consider a separate expansion for small 
and for large y. For y = 0(1) the previous expansion 
is valid and the probability distribution function in this 
region is given by 



pin ~ ^/ e -Su*-v\y\ ^ (33) 

with S given by equation (|23[) . This solution is sufficient 
to calculate the moment m provided that m < v. Indeed, 
we have 

( x m ) = [ e^y+^Pdydw = 5 - , X "> m /" . (33) 
J v A — m z 

Thus in this case, we obtain the deterministic scaling 
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B. 



The case v < n 



m/n for m < v and n < v. 



(34) 



For m > v the integral in equation (|33j) diverges and 
the calculation for the moments fails because the ex- 
pansion does not capture the large y behavior of the 
pdf. To remedy this we need to calculate the large y 
behavior of the pdf P out . We then rescale variables 
to y = y' — ln(/i)/n and w = w' + ln(/i)/n (so that 
y' + w' = 0(1) > 0), and obtain 

d w ,{e n{w ' +y,) P out ) + ud y ,P out + d 2 ,P out = . (35) 

Making the change of variables P out = e - vy '- nw 'Q and 



e ny /n, we arrive at 
d T & = -vd r Q + d r rd r Q 



(36) 



which is an advection-diffusion equation with space- 
varying diffusivity. The advecting velocity v = v In is 
directed away from the r = boundary. The boundary 
conditions are O — > for r 2 + r 2 — >• oo and 8 = /(r) for 
r — > and r finite, where /(r) is determined by match- 
ing with the inner solution. This problem can be solved 
exactly and its solution is given by 



9(r,r) 



1 



r(*) J- 



f(ro) 



0+1 



e -r/(r-ro) dTo ^ 



To obtain the functional form of / we match at an inter- 
mediate value of y chosen to be a/n ln(/x) with < a < 1. 
We obtain that to first order in fj, 



/(e—')-// -S(w'-w ) 



thus 



P° ut (y',w') 



vn \r 



e nw exp 



n(e~ 



2n y r(£>) (e~ 



(37) 

for w' < wo and zero otherwise. The calculation of the 
higher moments m > v can then be performed as 



-f%) 



— oo J — OO 



+oo /•+oo 



e my+mw P in (w,y)dydw 



e my+mw P out (w,y)dydw 

In the limit — > and for m > v the main contribution 
comes from the second integral resulting in 



zAT(m)XQ J(m, n) 
2r(£)n' > 



-m+l 



/' 



(38) 



m n .m—v—X 



du. Thus the seal- 



where J(m, f) = J °°(u + 1) 
ing exponent /3 m of the m-th moment is 



Pm = v/n for m > v and v > n. 



(39) 



When v < n the previous expansion fails due to the 
divergence of the denominator in cq. (|T9"|) . Figure H] 
shows the location in phase space of 10 3 trajectories for 
different values of [i (different colors) obtained by the 
numerical integration of the Langevin equations. It can 
be seen that as \x becomes smaller the distribution moves 
to smaller values of w but retains its width, unlike the 
case of the steep potential. We thus need a different 
expansion. Making the substitution w = w' + — ln(fx), 
the Fokkcr-Planck equation becomes 



CoP = fid' w [(l-e n ^' +y ^)P] 



(40) 



Since the derivative with respect to w' is multiplied by 
the small parameter, we write 

P = exp(^S(w'^j R(w',y), 

and expand R(w' , y) as R = Rq(w', y) + ^R\{w' , y) + 

At lowest order we obtain 

£ a R = o-( e "(«'+w) - 1)Rq + C R = 0. (41) 

where a = dS(w')/dw' . This equation can be solved 
exactly for positive and negative y. The two solutions 
are then matched at y = that selects the value of a 



2uI K Ue'™' /2 1 K K \\e nw 'l 



(42) 



where A 2 



-4a/n 2 



A 2 and L,- and K K 



are modified Bcsscl functions of order k. 

We can find an approximate solution of Eq. (|42|) to 
obtain a for w' — > — oo. To proceed we use the follow- 
ing relation for Bessel functions I a [z] K a [z] ~ 1/ (2a) + 
C\z la + C2Z 2 + . . . , valid for z — > and where C\ and C2 
are two constants. In this limit and for v < n, we obtain 
A oc e ™»V(20-»')), The asymptotic behavior of S is then 
of the form 



g n ( n ~ v ) nw'v/{n-v) 

Av 



(43) 



We observe that the exponential term e s ( w ^ M acts as 
a cut-off that selects very negative values of w' in the 
small-// limit. This limit will turn out to be useful when 
we calculate the moments. For the time being we proceed 
with our expansion without considering the w' — > — 00 
limit. 

The solution for Rq is Rq = A(w')e~ F ^ 2 g(w' ,y) where 
g is defined by 

K K \ e ^'+y)/2 j k Ae nu, '/2 ( y > ) 

(44) 

I K Xe n i w '+v)/ 2 k k \e nw '/ 2 (y < 0) . 

The solution Rq decays exponentially Rq ~ e uy for y 
—00 while for positive y the exponential decay (Rq ~ 
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e~ vy ) is followed by a super-exponential cut-off (R ~ 
cxp[-e ny ]) for y > -w'/(l-v/n). The amplitude A(tu') 
is obtained by a solvability condition on the equation at 
next order 



C a Ri = d w '[(l - e 



n(y+w') 



)Ro] 



(45) 



To obtain the solvability condition we need to multiply 
and integrate Eq. 05] by e F g(w',y) that is an element 
of the kernel of the adjoint operator of C a and thus the 
left hand side integrates to zero. The resulting solvability 
condition after some transformations reads 



d w , log ( AJ\e™'l n -Zo| ) = - T \ , 



2{e™'I n - 2b) 



where I q {w') is the integral 



e"«g 2 (w',y)dy. 



We can find the asymptotic behavior of Iq and I„ in 
the limit w' — > oo that leads to 1q ~ n 2 /2v 3 and 
I„ ~ (n/2v 2 )e~ nw . Inserting this result into Eq. |4T>1 
we obtain that 



A(w') ~ cxp 



2(n - v) 



w 



for «/ — > — oo. Already at this point we can observe by 
balancing the behavior of A with the expression for S in 
eq. HH that the most probable w scales like w <x ^ ln(/i) 
that is different from the scaling w cx — ln(/x) observed in 
section IV A. A comparison of the asymptotic result for 
the pdf with the results of numerical simulations can be 
seen in the lower panel of figure 2] 
We can now calculate the moments 



where we have introduced 



(47) 



M„ 



R {w' ,y)e qw ' +qy e s t w ">^dw'dy . (48) 



The small fi behavior of M q is obtained by keeping in 
mind that we can restrict to very negative values for w' 
which simplifies the expression of the Bessel functions. 
For m > v, the major contribution of the integral M m 
comes from the large y and large w part of the pdf that 
scale like y ^ —w ^ — l/fln(/i). For m < v the major 
contribution comes from the small y ~ and large — w ~ 
-in(/i). Careful evaluation of the integrals in the \i — > 
limit then leads to 



f3 m = min 



m 
v 



for v < n. 



(49) 



C. General expression for j3 m 

We have obtained several analytical expressions for the 
critical exponents that depend on the values of m, n and 
v. These expressions given in Eqs. (|34I39I49|) can be 
written in a compact form 



min [to, v\ 



(50) 



that divides the (y, m)-parameter space in four distinct 
regions with different scaling behaviors, as displayed in 
fig. [5] and fig. O Depending on the statistical measure 



m 
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FIG. 5. Phase diagram displaying the four regions where 
different scalings of the moments j3 m are observed. 



examined (degree of the moment to) different transitions 
can be observed. Increasing v we observe that the system 
transitions from on-off behavior /3 m = 1 (when v is small 
and the noise is dominated by the ^-correlated compo- 
nent) to the anomalous scaling (3 m — rajv (for to < n) 
or j3 m = v/n (for m > n) and finally to the deterministic 
one (when v is large and the distribution of y is narrow). 
We refer to the two intermediate scalings (for v between 
to and n) as anomalous because they do not follow nei- 
ther the mean field nor the on-off prediction. We note 
that the exponents are continuous functions of the pa- 
rameter v but not analytic. Therefore, they cannot be 
captured MS cl single Taylor series valid over the whole 
parameter space. 

For moments of small degree, the expression for the 
anomalous exponent mjv does not involve the nonlin- 
carity n while for moments of large degree, the exponent 
v/n does not depend on the considered moment. 

For fixed value of v and of the nonlincarity n, the scal- 
ing of the moments with m contains two regimes. A lin- 
ear behavior for small m and a plateau for large m. The 
value at the plateau depends on the width of the noise: 
corresponding to on-off ((3 = 1) for wide noise (small v), 
and a different value v/n for a narrow noise (large v). 
We point out that the moments do not depend linearly 
on to, i.e. that the solutions display multiscaling. This 
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traces back to the non-trivial expressions found for the 
p.d.f. obtained in section V. 

In figure |5] we display the first 4 exponents measured 
from numerical simulations. To obtain these exponents 



2.5 F 




0.5 
0.0 







O <* O 



V 



FIG. 6. Critical exponents of the first 4 moments (Pi dia- 
monds, 02 triangles, /?3 squares, 04 circles.) as a function 
of v for n = 2. The results were obtained from numerical 
integration and are compared with the analytical predictions. 



the Langcvin equations were solved for /x in the range 
1CP 6 < fi < 10" 3 and for a duration long enough for the 
4th moment to be converged. The exponents were then 
calculated by a linear fit. Simulations with (i as small as 
10 -7 for which only the first two moments were converged 
were also performed to verify that transient behavior as 
the one observed in section IV. C is not present. The 
measured exponents arc in agreement with the results 
obtained analytically. 



VI. A HEURISTIC DETERMINATION OF THE 
CRITICAL EXPONENTS 

In this section, we try to explain the predicted criti- 
cal exponents by giving a physical interpretation of the 
anomalous behavior. The evolution of z = log(a;) satisfies 
z = fi — x n + y. We note that close to criticality (fi <C 1) 
the amplitude of x remains small most of the time and 
thus the noise y is the dominant effect. Keeping only this 
effect, we have z — y that leads to the relation z = y + C, 
where the value of the integration constant C needs to be 
determined. Accordingly the marginal probability writes 
U z (z) = Il y (z+C) = jje~ v \ z - c \. This relation is only vi- 
olated at large x = 0(1) where the nonlinearities need to 
be taken into account and provide a large z cut-off. Tak- 
ing all these into account and returning to the x variable 
we can write the marginal probability B^ (x) as 



x- 1+u xZ. 



Il x (x) 



if 

if x c 
if x, 



< x < x c 

< X < Xj, 

,r < X 



(51) 



where x c = e c is still undetermined, x NL = 0(1) stands 
for the nonlinear cut off and N is determined by nor- 
malization. Figure [7] shows the marginal probability dis- 
tribution Tl x obtained from the numerical integration 
of the Langevin equations demonstrating that the two 
power laws model (|51[) gives a qualitative good descrip- 
tion of B^. Straightforward estimates of the moments 
(x n ) = J x n U x dx for this form of pdf result in 



„17l— V V 



if m < v 
if v < m . 



(52) 



x 
a 




FIG. 7. The marginal probability density function of £ , Tl(x), 
computed from the numerical solution of eq. [7] for v — 0.4 
and for three different /i (1CT 3 , 1(T 4 , 1(T 5 ). The position of 
the crossover between the power-laws moves to smaller values 
when fj, decreases. 



All that is left is to determine the value of x c - This 
can be obtained by balancing the averaged effect of the 
nonlinearity ((x n )) with the linear drift (/j,). As the es- 
timates for the moments indicate two cases need to be 
considered. 

For v > n, the average value of the nonlinearity is 
given by (x n ) cx a;™. Balancing with the drift leads to 

x c ~/i 1/n - 

For v < n, the tails of the pdf and the large x cut-off 
determine the averaged value of the nonlinearity given 
by (x n ) oc a;™"" ~ x". Balancing with the drift leads 
to x c oc ii x l v . 

Inserting the expressions of x c into eq. (JHU), we ob- 
tain all the behaviors predicted by eq. ((511)) . It is clear 
from these arguments that the tails of the y distribution 
control the presence of anomalous exponents. If y has a 
narrow distribution (F is a steep potential) x does not 
deviate far from the deterministic value and thus mean- 
field scaling is obtained. On the other hand, if y has a 
wide distribution it is the tails of the pdf, and the rare 
visits of a; to the nonlinear regime x ~ x NL that deter- 
mine the balance with the nonlinearity and the resulting 
scaling is anomalous. 
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VII. CONCLUSIONS 

Wc have introduced and studied the behavior of a fam- 
ily of zero-dimensional models in the vicinity of the in- 
stability threshold. The amplitude of the unstable mode 
in our models evolve in the presence of multiplicative 
fluctuations that makes it possible to have bifurcation at 
zero dimensions. Because this model is zero-dimensional 
(does not depend on space), it is among the simplest 
that can be considered. However, despite its simplicity 
the model exhibits nontrivial behavior. Depending on 
the control parameters the solutions display anomalous 
scaling close to the onset of instability: moments scale 
with the distance to onset as power-laws different than 
the ones predicted by mean-field theory. In addition, the 
system displays multi-scaling: the exponents are not sim- 
ply proportional to the moment degree. 

The values of the exponents were obtained by an exact 
calculation through perturbative expansions in the de- 
parture from criticality. This differs from what is usually 
obtained in equilibrium phase transitions where the ex- 
ponents are expressed as a series in the spatial dimension 
minus the critical dimension. In addition to the expan- 
sion we have presented heuristic arguments that allow 
to determine the critical exponents. This enables us to 
identify the basic ingredients for obtaining the anomalous 
behavior. First the model relies on a noise which spec- 
trum vanishes at zero frequency. We note that in the 
context of advection of a passive scalar in a turbulent 
flow, stochastic processes that have vanishing spectrum 
at zero frequency have been used to model anomalous 
diffusion |l4| . In the present model vanishing spectrum 
at zero frequency is achieved by considering the deriva- 
tive of a noise y that follows a random walk in a confining 
potential. If y has a narrow distribution, normal scaling 
is obtained. If y has a wide distribution, truly anomalous 
behavior takes place. 

All the analytical results were tested and verified by 
numerical simulations. Convergence of the estimated ex- 



ponents from the numerical results proved significantly 
difficult due to the appearance of intermediate power 
laws that are present in some limiting cases. These in- 
termediate power laws can contaminate the value of the 
exponents if sufficiently small values of [i are not investi- 
gated. Experimentally if such cross overs exist it may be 
difficult to distinguish them from true anomalous scaling, 
given the experimental limitations. 

As an example wc mention that in a recent experiment, 
the dynamo instability was observed in a turbulent flow 
of liquid sodium [l2j], [l3||. The first moment displays 
an exponent 0.78 located in-between 1/2 (as expected 
for cubic nonlinearities) and 1 (as expected for on-off in- 
tcrmittency) . The present model gives a possible expla- 
nation for the observed exponent since /3i can be larger 
than the mean-field 1 /n prediction and smaller than the 
on-off exponent 1. However the dynamo equations are 
more complicated than the model considered here: two 
fields (magnetic field and velocity field) are coupled and 
depend on space and time. Whether and when the dy- 
namo problem can be reduced to the model studied here 
remains an open question. 

Several additional investigations can be thought of. 
Other critical exponents can be defined and studied. The 
response to a constant field and the associated suscepti- 
bility are of particular interest. It is expected that simple 
scaling relations between the exponents exist (as in equi- 
librium phase transitions) and can be obtained exactly in 
the present context. Finally, taking space into account 
is also an attracting path. A possible attempt being to 
search for an expansion in the dimension since we have 
obtained a solution in the d = case. 
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